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The  linearized  equations  of  motion  for  a  ferry  vehicle  attempting 
rendezvous  with  an  orbital  satellite  (2-dimensional  Clohessy -Wiltshire 
equations)  are: 

x  -  2uy  =  u1 

(1)  2 

y  +  2ux  -  3u  y  =  u2  ' 

where  u.  and  u_  are  the  components  of  the  rocket  thrust  (the  mass  < 

1  1  S 

of  the  ferry  has  been  normalized  to  unity) .  The  amplitude  of  the  thrust  ;j 

is  | 

|u|  =  (u*  +  u§)*  ' 

I 

and  it  is  assumed  that  besides  the  initial  and  terminal  conditions  for  (l) 
the  additional  restraint 

(2)  |u(t)|<K  (0<t£T)  i 

is  imposed. 

We  wish  a  thrust  function  u  =  (u^,u2)  such  that,  starting  with 
an  initial  position  (xQ,y^)  and  velocity  (*Q>yg)  t  =  0,  we 

rendezvous  at  t  =  T  so  that 

(3)  x(T)  =  y(T)  =  x(T)  =  y(T)  =  0. 

We  ask  that  T  be  less  than  a  full  period  so  uT  <  2r.  If  it  is  assumed 
that  the  rate  of  fuel  consumption  is  proportional  to  the  amplitude  of 
the  thrust  then  the  total  fuel  consumption  will  be  (proportional  to) 

T  T 

f (u)  =  J  (u?  +  u^dt  =  5  |u(t)  jdt 
0  x  *  0 


(4) 
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and  our  problem  is  to  determine  a  thrust  function  satisfying  (2)  and 
(3)  (subject  to  (l))  which  minimizes  (4),  if  such  exists. 

We  let  M  be  the  set  of  transfer  functions  (thrust  functions  u 
for  which  (3)  holds  from  (l)  and  the  given  initial  conditions) .  We 
employ  the  norm 

T  T 

(5)  ||u|L  =  [5  |u(t)|2dt]^=  [I  (u?  +  u?)dt]^ 

0  0  1  * 

and  consider  the  Hilbert  space  X  of  u  =  (u^,u,>)  for  which 
is  finite.  In  the  topology  thus  defined,  M  will  be  a  closed  linear 
manifold  (see  below).  For  any  K,  let  be  the  subset  of  M  for 
which  (2)  holds;  we  assume  K  sufficiently  large  that  is  non-empty. 

Then 

Lemma  I:  The  functional  f(u)  actually  achieves  it8  minimum  on  Mj^. 

Proof:  We  observe  that  for  any  u,  u  e  X 

T  T  i 

|f(u)  -  f(u)  I  <  J  ||u |  -  |u||dt  <  $  |u  -  u|dt  <  T*  ||u  -u|L 
0  0 

(using  the  triangle  inequality  pointwise  and  then  the  Schwartz 

inequality) .  Thus  f  is  (strongly)  continuous  on  X.  We  also  observe 
that  if  a  +  p  =  1  with  a  >  0,  p  >  0  then 
T  T 

f(au  +  pu)  =5  |au  +  pu[dt  <  j  (a|u|  +  p|u|)dt  =  af(u)  +  pf(u) 

0  0 

(using,  again,  the  triangle  inequality  pointwise)  so  f  is  convex. 

Since  f  is  convex  and  (strongly)  continuous,  the  sets 

Sfl  =  (usff(u)  1  a} 
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are  convex  and  strongly  closed  in  X.  Hence  they  are  also  weakly  closed 
(c.f.,  [l]  p.  422)  and  (c.f.,[2]  p.  6l)  f  is  weakly  lower  semi -continuous. 

In  order  that  u  be  a  transfer  function  (u  e  M)  we  must  have 
x=y=x=y=  0  at  t  =  T  given  Xg,yQ,iQ>y0  and  the  differential 
equations  (l) .  The  solution  of  (l)  is  given  by 

t  - 

x(t)  =  a_(t)  +  $  fM(t  -  rfuCrjdz, 

1  0 

t  5 

x(t)  =  a_(t)  +  5  f  (t  -  t)u(t) d*c> 

(6)  2  0 

y(t)  =  a,(t)  +  5  f^(t  -  T)u(T)dT, 

J  0 

y(t)  =  a,(t)  +  5  f^(t  -  T)u(x)dT, 

4  0 

where 

a^t)  =  xQ  +  xQ(£  sin  ut  -  3t)  +  6yQ(ut  -  sin  ut)  +  j  yQ(l  -  cos  ut) , 

a2(t)  =  xQ(4  cos  ut  -  3)  +  6uyQ(l  -  cos  ut)  +  2yQ  sin  ut, 

a3(t)  =  2.  &q( cos  ut  -  1)  +  yQ(4  -  3cos  ut)  +  yQ  sin  ut, 

a^(t)  =  -  2  Xq  sin  ut  +  3wyQ  sin  ut  +  yQ  cos  ut, 

and 

f^(s)  =  (-  sin  us  -  3s,  -[l  -  cos  us]), 

2 

f  (s)  =  (4  cos  us  -  3  ,  2  sin  us) , 

f^(s)  =  (Jj[cos  us^  -  1] ,  sin  us) , 

f4(s)  =  (-  2  sin  us  ,  cos  us) . 
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Thus,  the  conditions  (3)  (which  define  M)  are  equivalent  to 
(3')  <u,F1>  =  (1  =  1,2, 3, 4) 

where  F1  =  /(t)  =  f^T  -  t)  and  A±  =  -  ai(T)  for  i  =  1,2, 3, 4. 

(We  have  used  <a,{5>  to  denote  the  scalar  product  in  X, 

T 

<<>,{!>  =  J  [a1(t)pi(t)  +  a2(t)p2(t)]dt). 

The  condition  (2)  (holding  for  almost  all  t  e  [0,T])  is  equivalent  to 


(2')  <u,g)  <  K  5  |g(t)  |dt 

0 

Certainly  (2)  implies  (2')  as 


(all  g  e  X) . 


T  T  T 

<u,g>  =  J  u(t)g(t)dt  <  5  |u(t)  |  |g(t)  |dt  <  K  j  |g(t)  |dt. 

0  0  0 

Conversely,  if  (2)  were  false  on  a  set  of  positive  measure  there  would  be 

a  K'  >  K  and  a  set  S  C  [0,T]  of  positive  measure  such  that 


|u(t)  |  >  K>  >  K 


If  one  now  took 
g(t)  = 

then  g  e  X  and 


(t  e  S) . 

t  e  S, 
t  /  S, 


T 

<u,g>  =  5  u(t)g(t)dt  =  J  u(t)  |^|-|  dt  =  5  ju(t)  |dt  >  K'  meas(S) 


1 

K  5  |g  jdt  =  K  5  IttI^  =  K  J  1  dt  =  K  meas(S)  <  K'  meas(S) 
n  <5  lul  5 


while 
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which  would  contradict  (2'). 


It  follows  that  u  e  Mj^  is  characterized  (except  on  subsets  of 
measure  zero)  by  the  conditions  (2')  and  (31)*  Then,  if  fu  }  is  a  net 
of  elements  of  weakly  convergent  to  u,  one  has 

=  limfl  <ua, f1)  =  lima  a1  =  a1  (i  =  1,2,3, 4) 

T  T 

<u,g>  =  lim  <u  ,g>  <  lim  K  J  jg|dt  =  K  5  |g|dt  (g  e  X) 

0  0 

so  u  is  also  in  and  is  weakly  closed*  Also  (2)  implies, 

for  u  e  M^., 

T  T 

||u|L  =  [5  |u|2dt]*  <  [J  K2dt]*  =  KT* 

*  0  0 

so  is  strongly  bounded  and,  hence,  weakly  compact  (c.f.,[  l]  p.  425) 


Thus,  in  the  weak  topology,  we  have  a  lower  semi -continuous  function 
f  on  a  compact  set  and  (c.f. ,[2  ]  p.  65)  f  must  achieve  its 

minimum  on  M^.  QED 


We  have  not  shown  that  the  u  e  which  minimizes  f  is  -unique 
but  will  call  any  transfer  function  which  minimizes  f  (subject  to  (2)) 
an  optimal  (or,  more  precisely,  K-optimal)  transfer  function.  Now  that 
we  know  such  functions  exist  it  is  of  interest  to  characterize  them. 
There  is  a  theorem,  the  so-called  ’'Maximum  Principle",  due  to  Pontryagin 
(c.f.,[  3 ] ,  pp.  66-68)  which  provides  a  characterization.  Introduce  the 
functions  (4q,  • .  •  ,’|> ^)  satisfying  the  linear  system 

=  "  RjU.u)]^, 


(7) 


♦0(T)  <  0 
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where  the  original  equation  (l)  has  been  written  as  a  first  order 
system  for  x  = 

(1')  \  (k  =  0,...,5) 

or,  in  detail, 

xQ  =  |u|  =  (u^  +  u2)^  =  Rq(x,u) 

^  =  x2  =  R1(x,u) 

(1")  x2  =  2ux^  +  u1  =  R2(x,u) 

*3  =  x4  =  R3(x,u) 

x^  =  -  2ux2  +  3u  Xj  +  u2  =  R^(x,u) 

*5  =  1  =  R5(x,u) 

(the  original  coordinate  functions,  x  and  y,  are  now  x^  and  x^)  and 

40  =  °  (*0(tj)  <  0) 

=  0 

(7')  *2  =  -  +1  +  2tJ»4 

=  -  3u^4 

*4  =  -  2^2  -  ^3 
♦  5  =  0. 

Construct  the  function 

H  =  H(+,x,u)  =  VkRk(i,u) 

(8)  =  +Q|u|  +  1'1x2  +  ^2(2ux^  +  ux)  +  ^3x^  +  $4(-  2ux2  +  3u2x3  +  u2)  +  ^ 

=  ^Q|u|  +  (cp,u)  +  ftjXg  +  2u^2x4  +  ^x4  +  +4(-  2ux2  +  3u2x3)  + 
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where  (p  =  (<P]_><P2)  * 

The  Maximum  Principle  says  that  a  necessary  condition  that  u  be 
K-optimal  (since  f  is  convex  and  is  a  convex  set  this  is  also 
sufficient)  is  that  there  exists  a  solution  of  (71)  such  that,  for  almost 
every  t  e  [0,T], 

(9)  H(^,x,u)  =  max  H(*,x,u). 

|u|<K 

We  observe  that  the  maximization  in  (9)  involves  only  =  ^qM  +  (<p,u)  . 
This  clearly  attains  its  maximum,  subject  to  (2) ,  if  u  is  equal  to  a 
multiple  of  <p/j<p|,  so  that  (cp,u)  =  |«p||u|.  Thus  the  function  HQ  =  w|u| 
where  w  =  4q  +  |?|>  is  to  be  maximized  with  u  =  a  qp/|«p|  and  0  <  |u|  <  K 
Hence: 

K  if  w  >  0, 

a.  = 

0  if  w  <  0. 

It  is  convenient  to  normalize  H  and  ^  by  taking  ^q(T)  =  -  1  (it  must 
be  negative)  so 

w(t)  =  -  1  +  M  =  -  1  +  (^  +  *|)* 

As  all  that  concerns  us  is  the  sign  of  w,  we  multiply  through  by  the 
positive  function  (l  +  |®j)  to  get 

(10)  W(t)  =  *2  +  ^4  "  1  =  ~  1 

and  i 

K  <p/M  W>0 

(11)  u(t)  = 

0  W  <  0. 
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Since  W  =  0  only  at  isolated  points  (unless  W  =  0,  which  is  impossible 
if  u  is  to  be  K-optimal)  we  do  not  oare  particularly  how  u  is  defined 
there* 

We  note  from  (11)  that,  for  a  K-optimal  transfer  function  u, 

One  always  wishes  either  to  be  firing  at  full  thrust  in  some  direction 
(W  >  0)  or  to  be  coasting  (W  <  0) . 

The  total  fuel  consumption  is,  then,  proportional  to 
T 

(12)  f(u)  =  5  |u|dt  =  K-t(u) 

0 

where  -c(u)  is  the  amount  of  time  during  which  the  rockets  are  to  be  firing 
for  a  K-optimal  transfer  function  u*  Since,  for  a  given  K,  any  K-optimal 
transfer  function  u  must  give  the  same  value  for  f(u) ,  t(u)  must  depend 
on  K  alone,  call  it  tk»  Thus,  for  any  K-optimal  transfer  function  u, 
the  total  amount  of  "firing  time"  is  and  the  total  amount  of  "coasting 

time"  is  (T  -  t^) ,  f(u)  =  Kr^. 

Clearly,  if  u  satisfies  (2)  for  any  K  it  also  satisfies  (2)  for 
any  larger  value  of  K.  Thus  C  M^,  for  K  <  K*  and 

(13)  Kt  =  min  f(u)  >  min  f(u)  =  K't  (K  <  K') 

ueMK  ueMK, 

so  that  not  only  is  xv  monotonically  decreasing  but 

(14)  tK  =  &(1/K) . 
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This  may  be  restated  as 

As  K  (the  maximum  available  thrust)  Increases  the  amount  of 
"coasting  time"  Increases  and  the  amount  of  "firing  time1’  goes 
to  zero  at  least  of  the  order  of  l/K .  One  expects,  "eventually", 
to  have  "infinitely"  powerful  thrusts  exerted  "instantaneously"  - 
l«e. ,  impulses. 

Of  interest  is  the  question  of  how  many  firing  intervals  there 
can  be,  which  is  equivalent  to  determining  a  bound  for  the  number  of 
times  the  rockets  can  "switch"  from  full  thrust  to  coasting  or  vice- 
versa.  From  the  characterization  (ll)  of  K-optimal  transfer  functions 
we  see  that  a  "switch"  corresponds  to  a  root  of  the  function  W  =  W(t) 
defined  in  (10)  from  the  system  (7') •  At  this  point  we  note  that  the 
general  solution  of  (71)  is  given  by 

*0=  cQ  (normalize  so  =  -  l) 


tj»  =  (c  .  -  3c.  t)  -  2c  „  sin  ut  +  2c,  cos  ut 
(15)  2  4  12  3 

^2  =  -  2u(c^  -  3Cjt)  +  3c2<j  sin  ut  -  3c^u  cos  ut 
♦  ^  -  (2/u)c^  -  Cg  cos  ut  -  sin  ut 


so  that  W  =  |cp  | 2  -  1  =  xtig  +  ^  -  1  has  the  form 

W(t)  =  3r2  sin2(ut  +  <{))  +  4r(c4  -  3c.jt)  sin(ut  +  0)  -  ^  c^y  cos(ut  +  (J) 
(16)  +  9c2t2  -  6c^c^t  +  (c2  +  c2  +  y2  -  l) 

It) 

l  c 

where  y  =  (c2  +  c2)^  and  $  =  -  tan  ^(~). 
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On  differentiating  (16)  we  obtain  (letting  0  =  ut  +  <J>,  for  convenience) 

_  2 
9  =  6y  sin  8  coe  6  -  8c^y  sin  8  +  4yu(c^  -  3o^t)  cos  6  +  lSc^t  -  6c^c^ 

and 

V  =  -  12y2u2  sin28  -  20c  1yj  cos  8  -  4yw2(c^  -  3c]Lt)  sin  8  +  (18c2  +  6y2u2)  . 

Letting  W  =  G  sin  8  and  writing  G  in  terms  of  Z  =  cos  8  gives 

G  =  -  12y2u2(1  -  ?2)^  -  20c1yu^(1  -  ?2)"^  -  4yu[(uc4  +  3^)  -  Jc-jOOS-1?] 

+  (18c2  +  6Y2u2)(1  -  ?2)-* 

and,  differentiating  with  respect  to  Z  (G1  =  dG/ d Z)  » 

G»  =  12y2w2C(1  -  ?2)“^  -  20c1yvj(1  -  ?2)" ^  -  20c]yu?2(1  -  ZZ)~3^2 

-  12yuc1(1  -  ?2)~*  +  (18c2  +  6  Y2u2)  ?(l  -  ?V3/2 

=  2y2u2(1  -  ZZ)~3^Z[W  -  ?2)  -  I6a(l  -  ZZ)  -  10a?2  +  3(1  +  3a2)?] 

=  -  2Y2u2(1  -  ?V3//2[6?3  -  6a£2  -  9(1  +  la2)?  +  16a] 

(where  we  have  introduced  a  =  c^/yu  to  simplify  the  expression).  Since 
G  =  G*?  =  -  u(l  -  ?2)^  G'  this  gives 

(17)  (1  -  ?2)G  =  12y2<4?3  -  a?2  -  |(1  +  a2)?  +  |  a.). 

For  any  value  of  the  parameter  a ,  the  cubic  polynomial  on  the  right 

of  (17) (and,  a  fortiori,  G)  will  vanish  for  at  most  three  distinct  real 
values  of  ?  (we  are  interested  only  in  those  roots  which  lie  between 
-  1  and  +  1,  however,  so  these  may  not  all  be  significant) .  Recalling 
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that  uT  <  2v  we  see  that  5  =  cos(ut  +  $)  can  take  any  value  at  most 
twice  for  0  <  t  <  T  and  there  is  an  upper  bound  of  six  on  the  number 
of  roots  of  &  in  [0,T].  Since  there  must  be  a  root  of  d  between  any 
two  roots  of  G,  it  follows  that  G  can  have  at  most  seven  roots  in 
[0,T].  Now  V  =  G  sin  ©  can  vanish  only  when  either  of  its  factors  does 
and  sin  S  vanishes  at  most  twice  in  [0,T]  so  V  can  have  a  maximum 
of  seven  plus  two  or  nine  roots  in  [0,T] .  As  above,  this  gives  a  bound 
of  ten  on  the  roots  of  W  and  of  eleven  on  the  possible  roots  of  W  in 
[0,T].  Each  "switch"  from  "full  thrust"  to  "coast"  (or  vice-versa)  of 
the  rocket  will  occur,  by  (ll),  only  at  a  root  of  W  in  [0,T].  Thus, 
we  have  shown  that 


There  are  at  most  6  firing  intervals  during  the  trip  for  a 
ferry  making  a  rendezvous  under  an  optimal  thrust  program. 


These  results  can  be  extended  somewhat.  Suppose  one  were  to  consider, 
instead  of  the  linearized  equations  (l),  the  equations  of  motions 


•  _  o  _  /  2  .  2\ 

X0  -  R0  -  (U1  +  U2' 

X 

0 

c+ 

11 

O 

— 

*1  =  h  =  x2 

xi(t0)  =  &1 

w = 

x2  =  H2  =  -  +  X3)  _3^Z  +  x"1  U1 

X2(V  =  a2 

x2(tl)  = 

(18)  x^  =  =  x^ 

X3(V  =  a3 

x3(t1)  = 

i4=R4=-x3(x?  +  X3)"3/2  +  X5_lu2 

x4(t0}  =  a4 

x4(tl}  = 

x$  =  R5  =  -  a(u^  +  Ug)* 

x5(t0)  =  mQ 

— 

*6  =  R6  =  1 

C+* 

II 

c+ 

0 

x6(ti)  = 

T. 


Here  (u^Ug)  is  to  be  chosen  so  as  to  minimize  xQ(tj)  subject  to 
meeting  the  initial  and  terminal  conditions  specified  in  (18)  and  the 
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same  condition  (2).  As  above,  x^t^  represents  total  fuel  consumption 
over  the  trip;  x^  and  x^  are  the  coordinates  and  x2  and  the 
velocity  components  of  the  ferry.  Also,  x^  is  the  mass  of  the  ferry 
(the  decrease  in  mass  of  which  with  fuel  consumption  is  now  being  con¬ 
sidered)  ;  x^  is  introduced  to  avoid  an  explicit  condition  on  t^. 

The  Maximum  Principle  now,  as  with  (7)  above,  introduces  new  functions 
{iji0, . . .  satisfying  the  system  of  equations 


40  =  ° 

\  =  (x2  -  2x^)(x^  +  X32)-5/2  *2  - 


*2  =  '  *1 

(19)  ^3  =  "  3xix3(xi  +  x^  ^  *2  +  ^*1  "  2x3^*1 

V  “  *3 

^5  =  “  x5  u1^2  “  x5  u2^*4 
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and  this  will  be  maximized  by  taking 


K<p/ 1  <p  | 

(22)  u  = 

0 


w  >  0 
w  <  0 


where  we  now  have 


(23)  w  =  |«p  |  -  m(l  +  atli^) 

with  =  (i|)„,tl>.)  as  before  and  m  =  x, .  Note  that  m  =  xc  and  a 
k  '  4  ?  5 

are  always  positive  and,  as 

*5  =  -  xj2  u  •  «p  =  -  x“2  |u !  |flp  |  <0,  +5(t1)  =  0, 

we  also  have  ^  always  non-negative  for  t  e  [tg,t^] .  As  before, we  can 
therefore  multiply  w  by  the  (always  positive)  function  [|«p|  +  m(l  +  a^)  ] 
to  get 

(24)  W  =  +  +4  -  “2U  +  a^)2 

K<p/|<p|  W>0, 

(22')  u(t)  = 

0  W  <  0. 

It  follows  that,  for  the  equations  of  motion  (IB),  the  K-optimal 
transfer  functions  (we  have  not  yet,  of  course,  shown  that  such  exist) 
are  characterized  in  exactly  the  same  way  as  for  the  linearized  equations 
( l) .  That  is 


For  a  K-optimal  transfer  function  (for  (18))  one  is  always  either 
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Unfortunately,  it  is  no  longer  possible  in  the  non-linear  case 
to  find  the  general  solution  of  the  auxiliary  system  (19)  independent 
of  x  and  u  nor  can  we  write  down,  as  with  (16),  a  simple  expression 
for  the  manifold  of  possible  W.  Thus,  the  method  used  above  to  demon¬ 
strate  a  bound  (independent  of  K)  on  the  number  of  firing  intervals 
cannot  be  employed  here.  We  have  no  proof  that  for  a  given  K  the 
determining  function  W  will  have  only  a  finite  number  of  roots,  nuch 
less  an  explicit  bound  independent  of  K. 

The  existence  theorem  given  above  for  the  case  of  the  linearized 
equations  (l)  no  longer  applies  as  the  manifold  M  cannot  be  charac¬ 
terized  as  in  (3')  -  the  proof  is  valid  modulo  a  lemma  to  the  effect 
that  M  is  weakly  closed.  It  is  possible,  however,  to  prove  existence 
by  the  use  of  a  general  existence  theorem  due  to  Roxin  [4]  which  for  our 
purposes  may  be  stated  as  follows: 

Theorem:  Let  x  =  R(x,u)  where  x,R  e  E  and  u  e  U  C  E  for  each  t. 

-  m  n 

Let  S-  =  [r(x,u):  u  e  Uj  be  convex  for  each  x  e  Em.  Assume  that  R  is 
defined  on  Em  x  U  so  as  to  be  continuous  in  u  and  Lipschitz  continuous 
in  x  (uniformly  for  u  e  U) .  Specify  initial  conditions  for  x(tg)  and 
terminal  conditions  x(t^)  e  E  where  E  is  compact.  Then  there  exists 
an  optimal  control,  i.e.,  a  measurable  function  u(t)  such  that  the  condi 
tions  on  x  are  satisfied  and  Xg(t^)  is  a  minimum. 

The  right  hand  side  R(  x,u)  in  (18)  is  not  uniformly  Lipschitz 
continuous , as  it  has  a  singularity  at  x^  =  x^  =  0  and  at  x^  =  0.  This 
is  easy  to  fix  up,  however,  -  if,  as  we  assume,  the  capsule  is  already 
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I 

j 


in  an  approximate  orbit  it  is  clear  that  the  trajectory  determined  by 
any  optimal  control  will  be  bounded  away  from  the  singularity  of  the  field 
(the  center  of  the  earth  -  actually,  of  course,  it  must  be  bounded  away 
from  the  earth's  surface)  and  x^  is,  of  course,  always  greater  than  the 
payload  so  that  R(x,u)  can  be  re-defined  within  that  bound  arbitrarily 
so  as  to  make  R  Lipschitz  continuous  in  x  without  affecting  the 
physically  significant  problem. 


That  S-  is  not  convex  for  each  x  (it  turns  out  to  be  part  of  the 
surface  of  a  cone)  is  a  more  difficult  matter  which  can,  fortunately,  be 
gotten  around  by  a  trick.  Consider  the  new  set  of  equations  (the  boundary 
conditions  are  the  same  as  in  (18)  except  that  we  place  some  a  priori 
bound  on  Xg(t^)  to  make  the  terminal  set  E  compact) 


(18') 


*1  =  *1  =  X2 

k2  =  R2=  ~  xl(xl  +  x3)_3/2  + 

*3  _  R3  _  x4 

*4  =  R4  =  -  x3^xi  +  *3> "3/2  + 

*5  =  R5  = 
x6  =  R6  =  1 


u 

u 

u 


1 

2 

0 


where  the  new  control  u'  =  (u^u^Ug)  lies  in  the  set  U'  C  E^  defined 
by 

(25)  0  <  (u*  +  u*)^  <  uQ  <  K. 


Modify  R',  of  course,  as  described  above  so  as  to  make  it  Lipschitz 
continuous  even  for  small  x^x^jX,..  Then  this  problem  satisfies  all  the 


i 
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conditions  of  Roxin's  existence  theorem  and  an  optimal  oontrol  u*(t) 
must  exist. 

If  we  could  show  that  necessarily  u^  =  (u£  +  u!p  *  for  this 
optimal  control  then  the  solution  to  this  problem  would  also  give  a 
K-optimal  control  for  the  problem  posed  by  (18) .  We  make  use  of  the 
Pontryagin  Maximum  Principle  again  and  note  that  (letting  il>  be  the 
solution  to  the  auxiliary  system  corresponding  to  (18')  -  which  happens 
to  be  the  system  (19)  again)  the  function  H  =  H(6,x,u')  is  now  given 
by 

(21')  H  =  [-  uQ  +  x“1(tl)2-u1  +  i!)^u2)  -  cnj^Ug] 

+  +  1>2X^)  -  (x^  +  *3)  _^^('l>2xi  +  +  ^6 

and  this  is  clearly  maximized,  for  u1  e  U,  by  letting 

?  ?  1  K  w  >  0 

(u^uj  =  u  <p/i<?|,  u_  =  (uf  +  uf)»  = 

1  *  U  u  x  0  w  <  0 

30  that  Uq  can  be  replaced  in  (18')  by  |u |  and  (18 ' )  reduces  to  (18). 

Thus  u  =  (u-^Ug)  provides  a  K-optimal  control  for  (18)  and  existence 

is  demonstrated. 
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